Time varying gauge couplings and di-nucleon states 



Malcolm Fairbairn 
Service de Physique Theorique, CP225, 
Universite Libre de Bruxelles, B-1050 Brussels, Belgium 

In this talk we discuss the effect upon the di-proton, di-neutron and deuteron of a time 
variation of the QCD gauge coupling. We describe how a time evolution of the size of extra 
dimensions can give rise to such a variation and show how this is related to the recent results 
on possible time variation of the electromagnetic fine structure constant. Based on work in 
collaboration with Tom Dent. 



1 Introduction 

Webb et al have observed the absorption of light from quasars by intervening intergalactic 
clouds of matter at high redshiftsS. By studying the relative position of absorption lines in 
these clouds, they have found that there appears to be a variation in the electromagnetic fine 
structure constant a em at high redshift. 

Such a variation is naturally accommodated in higher dimensional theories such as string 
or M-theory since in these models the value of the gauge couplings are set by the size of the 
compactified directions. Cosmological evolution of the higher dimensions would then lead to 
time variation of the gauge coupling of the grand unified theory. One would expect this to lead 
to a variation of all of the gauge couplings in the standard model including the QCD coupling. 

The renormalisation group running of the QCD coupling constant 0:3 means that the accom- 
panying fractional variation in the characteristic energy scale of hadronic physics Aqcd would 
be much larger than that for a em . It is natural therefore to investigate as to whether or not the 
observed variation is consistent with the assumptions of gauge coupling unification and physics 
in the hadronic sector. 

After briefly discussing the results, we talk more about the relationship between the size of 
higher dimensions and the gauge couplings in the 4D theory in one particular model. Next, we 
describe how variations in a em lead to variations in Aqcd and the implications of this for the 
binding of the two nucleon states, the deuteron, di-neutron and di-proton. We then show that 
this leads to constraints on the variation of a em from the variation of the binding energy of the 
deuteron during nucleosynthesis. 



2 The observations 

As we have said, Webb et aS have measured a em at high redshifts by observing the absorption 
lines in the light from high redshift quasars. This absorption occurs in clouds of intervening 
intergalactic dust. When trying to understand the origin of a time variation in a em , it is useful 
to know a tiny part of the atomic physics that allow one to probe such a variation. 



Upon solution of Schro dinger's equation in atomic physics, the nth electron energy level of 
a given atom is proportional to 

me Z 

where m is the electron mass, e is it's charge and Z is the number of protons in the nucleus. 
The relativistic correction to this energy level has the following proportionality i 
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Comparison of the previous two equations is the key to understanding why observation of spec- 
tral lines in high redshift objects is an extremely sensitive probe of the variation of the gauge 
couplings. Typically, if one was trying to determine the energy of a particular spectral line 
in a high redshift object, there would be many systematic errors due to uncertainties in one's 
knowledge of the redshift or proper motion of the object. However, these systematic errors fall 
out of the analysis if one deduces a em by comparing a number of lines within the same cloud. 
It is also interesting to note that such measurements are not sensitive to a time variation in the 
mass of the electron, since this would simply act like an apparent change in the redshift of the 
system. 

It is then possible to compare the observed lines to those predicted numerically to see if 
the best fit. is with a em as we know it today, or with a different value of a em . The analysis of 
Webb et a£l shows that for redshifts of approximately z ~ 1 — 3 the observed spectral lines are 
statistically significantly better fit by the numerical models if one assumes that a em was smaller 
by about 1 part in 10 5 at those times. 

It is therefore interesting to consider the effect of this variation assuming that the gauge 
fields are part of some larger GUT symmetry and that variation of one gauge coupling leads to 
variation of all three gauge couplings. 

3 GUT theories and higher dimensions 

The observed unification of the gauge couplings at 10 16 GeV in the supersymmetric (SUSY) 
version of the standard model stronglp suggests that both QCD and the electroweak theory 
are parts of a larger gauge symmetry broken at high energy. In this section we show how the 
coupling of the GUT gauge field is related to the size of extra dimensions in one particular 
model, namely Heterotic M-theory. We then report how this all relates to a em and the hadronic 
sector following the analysis of our papeB. 

3.1 Some (M)otivation 

We know that in order to have general relativity and some unified gauge theory with field 
strength F containing the standard model, the Lagrangian of the effective four dimensional 
theory obtained by compactification of some particular string/M-theory must look like 
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As we mentioned earlier, the values of the gauge and gravitational couplings oigut and G 
are set by the size of the extra dimensions. This arises because the volume of the higher 
dimensional space in which each particular field propagates re-scales the coupling of that field 
upon compactification. 

For example in Heterotic M-theory, the GUT gauge field comes from the Eg gauge field in 
the 11 D theory a. This gauge field exists only in the uncompactified directions and in the 6 



extra directions which are compactified upon a Calabi-Yau manifold. The coupling in 4D is 
given by 

2vr 2 / 3 K 4 / 3 
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where k 2 is the 11D Newton's constant and V is the volume of the Calabi-Yau. As usual, gravity 
propagates in all space-time directions including the 11th dimension and the effective Newton's 
constant in 4D is given by 

k 2 

G = mWp (5) 

where p is the size of the 11th (orbifolded) dimension. So it is possible to have variations in the 
4D strength of gravity if p evolves over time and in the strength of both gravitational and gauge 
interactions if V evolves over time. 

There have been various suggestions as to the mechanisms responsible for stabilising V 
and p at least in weakly coupled heterotic string theory. Perhaps the most well known are 
non-perturbative corrections to the scalar potential for these fields due to the condensation of 
the fermionic superpartners of gauge bosons oQ, but the problem is not completely solved and 
becomes more complicated in the strong coupling regime. 



3. 2 Other aspects of time varying extra dimensions 

In all compactifications of higher dimensional theories, if the sizes of the extra dimensions are 
allowed to possess dynamics they give rise to scalar fields in the 4D theory with cosmological 
equations of motion of the usual form 

4> + 3H<j) + m 2 <p = (6) 

where m is the second derivative of the potential in which the field is rolling and H = a/ a is 
the Hubble expansion factor. Variation of the size of the higher dimensions over time leads to 
variations in these scalars 4> and to the variation of our effective 4D gauge couplings. Now the 
observed variation of the fine structure constant is of order Aa em /a ern ~ 10~ 5 at a redshift of 
order z ~ 1@. Nucleosynthesis occurs at redshifts of order z ~ 10 8 — 10 10 and the variation in 
gauge couplings at that time is certainly constrained to at least Aa em /a em < 1. If this were 
not the case, it would not be possible to obtain a primordial ratio of Helium and Hydrogen that 
matches that observed in regions of the universe where no nuclear processing has occurred since 
nucleosynthesis™. Therefore if one is to take seriously dynamics of the form of equation @ 
it is either necessary to cook up a very original potential, or to assume that the slope of the 
potential is extremely small so that equation <j) is damped and the <ft is frozen at early times 
when 

m 2 <C H 2 . (7) 

The variation of the gauge coupling would then begin later on when the expansion drops to the 
point where H 2 ~ m 2 and the field <f> is no longer frozen. This leads to a fine tuning problem 
in explaining why m is so small Also one might expect a very small mass for such a field to 
lead to dangerous light modes 113 and/or have associated with it a particle leading to new long 
range forces 1111. Therefore a variation at such a late time, if proved to be correct, would create 
so many problems for our existing model that it would give us a great deal of information about 
not just unification but many other aspects of particle physics beyond the standard model. 



3.3 Aqcd = f(oi e m) i n GUT theories 

The relationship between Aqcd an d the gauge coupling in the grand unified theory is obtained 
by considering the relationship between Aqcd and the QCD gauge coupling at energies just 



above the masses of the light quarks. Then one can find the relation between at these low 
energies and at the energy scale of unification!! Mjj to obtain 



_2_ 

m c mbmt\ 27 / —6tt 



A OCD = Af * AT,? c 6 f exp f— (8) 

where is the mass of the i quark and we have assumed that all the SUSY particles have 
the same mass which is equal to the mass of SUSY breaking Msusy- Next, by running the 
electroweak gauge couplings a\ and a.2 down from Mjj to low energies we obtain the relation 
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the linearity being applicable for variations in a em less than a few parts in 10 3 . 

It is worth pointing out that at high energies, and therefore high temperatures in the early 
universe, one would expect the same renormalisation to change the value of the gauge couplings 
and hence a em . How do we know this is not responsible for the possible time variation of a em ? 
Well apart from the fact that the spectral lines correspond to a jump in electron energy and 
therefore are more like a tunnelling process than something like compton scattering, the variation 
goes the wrong way and one would expect no renormalisation of a em below the electron mass 
(1 MeV, z~ 10 10 ) anyway 



4 The effect of varying gauge couplings on di-nucleon states 

4-1 The internuclear pion potential 

If one ignores the small contribution due to the bare quark masses and the interquark electro- 
magnetic binding energy, the nucleon mass Mjv and the masses of mesons such as the sigma and 
omega meson scale with Aqqd 

M N cx m w oc m a oc Aq C d- (10) 

This is not true for the pion mass which arises as the Goldstone boson of spontaneously broken 
chiral symmetry. Goldstone bosons of exact symmetries are massless, but the fact that the 
underlying chiral symmetry of the QCD Lagrangian is already slightly broken by the presence 
of small but non zero quark masses leads to a pion mass which scales as (see e.g. El) 

1 
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where m q is the average mass of the light quarks. Here we have also assumed that the energy 
scales of confinement and chiral symmetry breaking coincide as suggested by simulations on the 
lattice The exchange of pions contributes to the internuclear potential as 

V(r) = - J —- (12) 
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where the strength of the force is given by 

f 2_( 9>4\ m) 

where g n is the pion nucleon coupling and is the nucleon mass@. 



4-2 The stability of the di-nucleons 



The fact that the di-proton and the di-neutron are not bound is essential for the universe to 
evolve the way it has done. This is because the existence of a stable di-proton would lead to 
a rapid channel for the burning of hydrogen in stars. If the di-proton was stable, it would not 
be possible for stars to slowly burn hydrogen on the main sequence for long periods of time so 
the fact that the di-proton is unbound is essential for our understanding of stellar evolution. 
Since the internuclear potentials in the di-proton and di-neutron systems are described by the 
same nuclear potential (|l2|), the only difference between the two binding criteria arises from 
electromagnetic effects which are very much smaller than the precision of the results in this 
paper. 



We can use the potential (12) to calculate the amount that one would have to change the 
ratio between m q and Aqcd in order to stabilise the di-neutron or di-proton system. We use a 
trial wavefunction of the form 



so that at small r the second derivative of the wavefunction depends mostly upon the potential 
whereas at large r it depends mostly upon the energy, denoted by the variational term b. Then 
we can write 




Figure 1: Stability of di-nucleon system. Dark region is bound, light region unbound. 
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Once we have evaluated v, we can see from figure |l] that in order for the di-nucleon system 
to remain unbound, Av/v < 2.6. This corresponds to a very large change in Aqcd and the 
corresponding change in a em is about -5%. A similar analysis for the deuteron leads to a 
criterion for un-binding of about a +2% increase in a em . 



4-3 Nucleosynthesis Constraints 



We have placed constraints on the amount by which it is possible to vary a em so that the 
corresponding variation in Aqqd is small enough to be consistent with the stability of the 
deuteron and the instability of the di-nucleon. Of course, without saturating these bounds, 
a variation in the ratio between the quark masses and Aqcd will lead to a change in the 
binding energy of the deuteron. This will have implications for the production of Helium during 
nucleosynthesis. 

When the universe is hotter than about 1 MeV, neutrons and protons are continuously 
interconverting via the weak interactions 

n + e + <-> p + v e 
n + u e <-> p + e~ 

n <-> p + e~ + P e . 



with the rate of these reactions being given approximately by 

^pe-^nu ~ G^T 5 (16) 

Now when the Hubble expansion H = 1.66y/g*GT 2 becomes larger than this rate the reactions 
freeze out. This occurs at a temperature 

Ql/6 
Lrp 

and the ratio of neutron to proton number density (n/p) which one finds at this temperature is 
frozen, only changing via the beta decay of the neutron. 

(n /p ) = e-( M "- M ?y T ° (18) 

Later at T# e ~ Bd eu t, the deuteron can survive without becoming photo-dissociated by the hot 
plasma and the following Helium forming reactions are able to take place 

H 2 (H 2 ,ra)He 3 (H 2 ,p)He 4 
H 2 (H 2 ,p)H 3 (H 2 ,n)He 4 

H 2 (H 2 , 7 )He 4 . 

Essentially all (99.99%) of the neutrons around at this temperature form helium via deuterium. 
The mass fraction of Helium is then given by 



Y = Xi= 2(„/p) 
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Now as the neutrons decay, the neutron proton ratio changes 



; x A = 2dd (i 9 ) 
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V ' ' (p /n ) + l-exp(-t/T) V ; 

where r is the neutron lifetime, so the time between weak freeze out and the temperature at 
which Helium can form is essential in determining the resulting mass fraction of Helium in 
the universe. Since this temperature is determined by the binding energy of the deuteron, the 
process is extremely sensitive to variations in the internuclear potential and consequently the 



gauge couplings. On finds that in order to agree with the observed primordial mass fraction of 
Helium the binding energy of the deuteron is constrained to lie within 



- 0.20 < - < +0.13 
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which corresponds to a variation in Aqcd of 
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If we calculate what variation in the gauge couplings can give rise to this variation in Aqcd we 
see that 

- 0.005 < < 0.005 (23) 

"3 ITeV, no SUSY 

for the variation in 03 at 1 TeV, assuming that is where the gauge couplings are set. If we 
consider the more conventional case with SUSY and gauge coupling unification at 10 16 GeV, we 
can find out exactly what this means in terms of the variation of cv em 
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(24) 



This shows us that the constraint on a em is very strong around nucleosynthesis, so if it is varying 
today, it cannot have varied much more rapidly in the early universe. 



5 Conclusions 



We have shown that a variation of a em implies a much larger variation in Aqcd if the gauge 
couplings are unified at very high energies. This leads to constraints upon the variation in a em 
due to the fact that we must maintain the stability of the deuteron and the instability of di- 
proton at all times since nucleosynthesis. These constraints are typically of the order for a few 
percent. We have also shown that it is possible to place new constraints upon the variation of 
a em at nucleosynthesis since quite a small variation in a em would give rise to a large variation 
in the primordial Helium abundance if one assumes gauge coupling unification. 
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